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Abstract. We study the valuation and hedging problem of European options in a mar- 
ket subject to liquidity shocks. Working within a Markovian regime-switching setting, 
we model illiquidity as the inability to trade. To isolate the impact of such liquidity con- 
straints, we focus on the case where the market is completely static in the illiquid regime. 
We then consider derivative pricing using either equivalent martingale measures or expo- 
nential indifference mechanisms. Our main results concern the analysis of the semi-linear 
coupled HJB equation satisfied by the indifference price, as well as its asymptotics when 
the probability of a liquidity shock is small. We then present several numerical studies of 
the liquidity risk premia obtained in our models leading to practical guidelines on how to 
adjust for liquidity risk in option valuation and hedging. 

Keywords: liquidity shock; indifference price; exponential utility maximization 

1. Introduction 

The traditional option pricing theory relies on a key assumption: that markets are 
always liquid and agents can trade whenever they wish. However, this is not generally true 
in real financial markets |16j . There is increasingly a common phenomenon of liquidity 
crises in even well-established securities markets, where liquidity dries out or trading is 
virtually impossible for a number of reasons, such as political turmoil, war, and financial 
crises. Two prominent examples include the terrorist attacks of 9/11 which closed all US 
security exchanges for four days, and the "flash crash" of May 6, 2010 [13] when numerous 
stocks experienced extreme price swings (up to 99% of value) coupled with minimal trading 
volume during the course of several hours. 

In recent years a number of approaches for dealing with market illiquidity within op- 
timal investment frameworks have been developed. Schwartz and Tebaldi [21] considered 
the optimal portfolio problem with permanent trading interruptions; Diesinger et al. [5] 
addressed terminal wealth maximization for CRRA utilities. Ludkovski and Min [17] and 
Gassiat et al. [9] analyzed the impact of illiquidity on optimal consumption strategies. In 
this paper we extend this analysis to the problem of pricing of derivative securities. To 
our knowledge, there is no previous analysis for effect of illiquidity on option prices. 

In this paper, we define illiquidity as the inability to trade in a timely way. While 
liquidity crises involve a host of phenomena that affect the market environment, includ- 
ing downward jumps in asset prices, modified asset dynamics and restrictions on trading 
strategies, here we concentrate on the last item. Namely, to isolate the impact of liquidity 
shocks, we consider the extreme scenario whereby (i) trading in assets is completely sus- 
pended during the shock and (ii) the underlying asset price remains frozen while the shock 
is ongoing. 

We adopt a regime-switching description of market liquidity which is modeled through 
a separate Markovian liquidity factor. This aspect of our model connects with the growing 
literature on Markov-modulated financial markets driven by the recognition that the eco- 
nomic environment is not stable. Some of the pertinent analyses include [H El [TTj [T5l [23] . 
The model of Blanchet et al. [3] effectively corresponded to the case where the liquidity 
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shock was permanent, while Siu et al. [23] and Leung [15] studied valuation of options 
using Esscher and minimal entropy martingale measures, respectively. 

The main rationale for our choices is to have a well-understood benchmark for compar- 
ing derivative valuations and to focus on the effect of trading suspension. In particular, 
thanks to the special form of liquidity shocks, we can directly compare to the classical 
frictionless Black-Scholes model, specifically linking to the options' time-decay. Indeed, a 
liquidity shock induces an instantaneous (in terms of business-time) jump in the option's 
Delta, making perfect hedging impossible and exposing investors to charm risk, namely the 
dependence of Delta on time-to-maturity. In reality, asset valuations continue to change 
(often in an unfavorable way to the investor) and some form of trading may still be possible 
(e.g., only a short-selling restriction). However, the resulting model is simply somewhere 
between the classical textbook description and the one considered herein; similarly, incor- 
porating other features, such as jumps in asset prices, only obscures the final conclusions. 

Potential liquidity shocks introduce a new non-traded source of risk and make the asset 
market incomplete. Derivative pricing in incomplete markets largely follows two main 
approaches. First, one may fix an equivalent martingale measure (EMM) Q and calculate 
option prices as expectations under Q. Some of the common candidates for EMMs include 
minimal martingale measure [7], minimal entropy martingale measure [U [19], variance- 
optimal martingale measure [22] . and empirical martingale measure [3]. Second, one may 
use a utility maximization criterion to obtain indifference prices that take into account the 
risk aversion of the investor [21 [TO] [TT| [T2l [T5l [20| [26] . Given the widely reported crash-o- 
phobia of liquidity crises, this is a useful nonlinear pricing rule complementing the classical 
expectations-based approach. In this paper we compare both approaches for determining 
the option's liquidity premia, with a special focus on exponential utility valuation, which 
leads to wealth-independent indifference prices and natural connections to the minimal 
entropy measure. 

Our main contributions concern the analysis of the Hamilton-Jacobi-Bellman (HJB) 
equations satisfied by the indifference prices, including explicit asymptotics in terms of the 
two main model parameters: level of risk-aversion and probability of liquidity shock. We 
also present a detailed comparative analysis between the model prices and the classical 
Black-Scholes prices with the aim of providing rules of thumb for incorporating liquidity 
shock risk into the valuations. For this purpose, we introduce the new concepts of implied 
and adjusted time to maturity that allow translation of liquidity risk premia into the Black- 
Scholes language. 

The rest of paper is organized as follows: in Section [2] we setup the probabilistic model. 
In Section [3] we summarize the EMM approach to pricing under liquidity shocks; Section 0] 
then describes the utility indifference approach to derivative valuation. Section [5] discusses 
asymptotic analysis of the utility-based pricing mechanism to gain further insight into our 
formulas; in Section [6] we present a handy interpretation of the liquidity premium through 
implied time-to-maturity. All these results are illustrated in Section [7] with numerical 
experiments. Most proofs are delegated to the Appendix. 

2. The Market Model 

In this section, we describe the market model in the presence of potential liquidity 
shocks. Fix a complete probability space (0, J-, P), where P is the real-world probability. 
We consider a finite investment time horizon T < oo, which is chosen to coincide with the 
expiration date of all securities in our model, and a financial market which consists of a 
stock and a risk-free asset (cash), and is subject to the liquidity shocks. We use (X t ,irt, St) 
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defined on 1 x K x R + to denote the total wealth, stock holdings as proportion of total 
wealth, and the stock price at time t £ [0, T], respectively. 

2.1. Market Dynamics. We assume the market has two states, liquid (0) and illiquid 
(1). We use a continuous-time Markov chain (Mt) to represent the changing liquidity of 
the financial market with state space E = {0, 1}. The Markov chain (Mt) is a proxy for 
the fluctuating market liquidity and modulates asset dynamics. Its infinitesimal generator 
is 

A=(~ U01 m ), (2.1) 

where vq\ and v\q are for simplicity constants. Associated with (Mt) are two counting pro- 
cesses (iVoi(t)) and (N\o(t)) with intensity rates of ^oil{A/t=o} an d ^lol-nw^l}; respectively 
counting the transitions — > 1 and 1 — > 0. 

We assume that the market value of the stock follows a Markov-modulated Geometric 
Brownian motion (GBM) model. Without loss of generality, we assume a zero interest 
rate which is equivalent to working with discounted price processes. In the liquid state 
(Mt = 0) the market dynamics follow the classical Black-Scholes model, and the stock can 
be traded continuously and frictionlessly. More precisely, we take 

dS t = fi S t dt + a S t dWt, (2.2a) 
dX t = fiomXtdt + a 7r t X t dWt, (2.2b) 

where //o and <ro denote the stock drift and volatility and (Wt) is a standard one-dimensional 
Brownian motion under P, which is assumed to be independent of the Markov chain (Mt). 
In the illiquid state (Mt = 1) the market is static and trading in stock is not permitted: 

dS t = dX t = 0. 

As can be seen, this market model has two sources of uncertainty, one generated by (Wt) 
and the other by (Mt). The liquidity constraint is captured by our set A of admissible 
trading strategies. Namely, letting Qt '■= o~(S s ,M s : s < t), a strategy (irt) € A is admis- 
sible if it is self-financing, ^-progressively measurable, satisfies the integrability condition 
^Uo d£\ < oo, and the constraint d%t = on {t : Mt = 1}. 

We focus on pricing European contingent claims with maturity date T and payoff 1i(St). 
For simplicity, we assume that the terminal payoff does not depend on Mt; it would be 
straightforward to incorporate liquidity penalties for option exercise. For later use we recall 
the complete Black-Scholes market counterpart which has the unique martingale measure 
P with 

dS t = a S t dWt, (2.3) 

and where contingent claim prices are given by the classical no-arbitrage formula (note 
that we use the time-to-maturity parametrization) 

P B s(T,S)±t[h(S T )\S Q = S], (2.4) 

where E = E p . 

Remark 2.1. More generally (see [5l[T7]) the market need not be static during the illiquid 
periods. Thus, one may postulate modified dynamics for the stock during the liquidity 
shock and also include instantaneous price drops when the liquidity regime changes. For 
instance, assuming that when Mt = 1 the stock follows a GBM with drift fj,\ and volatility 
o"i and that it experiences a fixed drop of L% at the beginning of each liquidity shock, the 
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analogue of flgSJ) would be S t = (1 -L)S t -, X t = (l-Lnt-)X t -, and vr t = vr t _(l - L)/(l - 
TTt-L) whenever Mt — Mt- = 1 and 

dS t = ii x S t dt + atS t dW t , (2.5a) 

dx t = nmx t dt + cr^tXt dw t , (2.5b) 

d-K t = {/il7T t (l - 7T t ) + ^ 2 t {2TT t - 1)} <tt + <Tl7T t (l - 7T t ) dW t . (2.5c) 

However, the resulting model is significantly more complex due to the fact that nt is chang- 
ing during the illiquid regime and hence remains part of the system state. In other words, 
the state variables when Mt = 1 are (irt,Xt,St) introducing an extra state dimension. We 
further discuss the resulting pricing equations in Remark 14.81 

Notations. We use Kj Si ,i = 0,1, to denote expectation taken under the measure P with 
starting values St = S, Mt = i. We also let 

= { °° = °' (2.6) 




D = ^ , (2.7) 

\o oj 

so that j3's are the Sharpe ratios in each state (note: the zero risk- free rate). 

3. Option Pricing via Martingale Measures 

In this section, we summarize the option valuation problem in the context of the EMM 
methodology. We mainly follow Leung [15] who considered a related regime-switching 
market. We begin by characterizing the set £ of all equivalent local martingale measures 
(EMMs) in this liquidity switching market. A probability measure Q in £, equivalent to 
P, should be regarded as a risk neutral measure with respect to both asset price and the 
liquidity risk. Using Girsanov's Theorem for the Wiener process (Wt), and the Markov 
Chain (Mt) we obtain (cf. Theorem 7 in Leung [15]): 

Proposition 3.1. The set £ of all possible EMMs is given by 

£ = |Q a : 3{o!t(i, j)}i^j,i,jeE > bounded and ^-adapted s.t. 

1p~L = ^ X ef EP Ko(*)Cf(*)] = l,Vt G [0,T]}, where 

£o(i) = exp (~ p 2 (M s )ds - J* f3(M s )dw)j , 

ef(t) = expf- f T (A s (M s ,M s )-A s (M s ,M s ))ds)- TT a s (M s .,M s ), (3.1) 

^ io / 0<s<f. 

M s -jtM s 

•n, ac -\ J a t(hj)A(i,j) if i^j, . 
with MhJ) = [_< LkfU A t (i,k) if i=j. and A given m (EU). 

In Proposition 13. 1\ the exponential martingale (£o(t)) is the Girsanov transform that 
makes (St) a Q Q -martingale, while (£f(t)) is the exponential martingale that transforms 
the generator of (Mt). To preserve the Markovian property of (Mt) under the new measure 
Q a , we henceforth require that all (at(i,j)) be Markovian. Then, under Q a , the generator 
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matrix of (M t ) is A = [A t (i, j)]ij & E- The collection {a t (i, j)}i^j can be considered as the 
risk premium factors for the liquidity switching risk |15j . 

Using the law of iterated expectations by expressing the price as a discounted risk neutral 
expectation over all possible paths of the liquidity state (Mt) we obtain 

Proposition 3.2. Given a measure Q a € £, the option price p a (t,S,i) is given by 

p a (t,S,i) = Ef tSA [h(S T )] = e -f t T ^(u)dup Bs{T _ tjfif) 

+ J h,l-i{r)e- ^ '^-^ du n,s,i\p a (r,S T , 1 - i)]dr, (3.2) 

where u id (t) = A t (i,j) = a t (i, j)A(i, j), i ^ j. 
A formal proof of Proposition 13.21 is given in [AJ 

In the special case where a = 1 is a unity transformation, we obtain the Minimal 
Martingale Measure (MMM) Q° (with ^| = £o(*)) which corresponds to a zero 

risk premium associated with the liquidity risk. By Girsanov's theorem, under Q°, (St) 
follows 

dS t = l{ A / t= o}Cr St dW t , (3.3) 

where (Wt) is a standard Brownian motion under Q° and the Q°-generator of (Mt ) remains 
A. We use p MM (t,S) to denote the price under MMM conditional on Mt = 0. Another 
special case is furnished by the minimal entropy martingale measure (MEMM), Q E . 

Definition 3.3. The minimal entropy martingale measure, Q E , minimizes the relative 
entropy with respect to P over the set of EMMs £ : 

Q E = argmm Qe£ H T (Q\F), where (3.4) 



H T (Q\P) = { ^ r B dP 1 V1 J ' ^ < ' (3.5) 

+oo, otherwise . 

Definition 3.4. Consider the ODE system F'(t) = (D — A)F(t), with F(t) = (F 1 (t),F 2 (t)) 
and F(T) = [1, 1] T . The solutions Fo(t), F±(t) are explicitly given by 

F (t) = Cl e Xlt + c 2 e X2t , (3.6) 

Fi(t) = — {ci(dv + z/ i - Ai)e Al * + c 2 (d + v m - \ 2 )e x A , (3.7) 
Urn { > 



2al 



where d = D(l,l) = ^, (3.8) 



d + foi + ^io ± \J (do + z^oi + ^io) 2 - 4d ^io 
Ai,2 g , 

Cl = ^0 e -A 1 T > and C2 = Al^^T 

A 2 — Ai Ai — A2 

Remark 3.5. Let F 2 (i) = e~ d °^ T ~ t \ Through simple algebra, we have < A 2 < do < Ai, 
which leads to iq(i) > F (t) > F 2 (i), Vt € [0,T]. 

Applying Theorem 9 in Leung [15] we immediately obtain the following characterization 
of Q E . 
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Proposition 3.6. The MEMM is Q E = Q qE , where the minimal entropy risk factor a E 
is given by a E (i,j) = £ {0, 1}. Thus, under Q E the Markov chain (M t ) becomes 

a time-inhomogeneous Markov chain with transition intensities 

h,i~i(t) = iv-i%^, 2 = 0,1- (3.9) 

Proposition 13.61 presents an explicit form of the MEMM. By Remark 13. 5| under Q E , 
(Mt) is more likely to be in state 1, z>oi > ^oi^io < ^io- That is, under the MEMM <Q> E , 
the market is more vulnerable to a lengthier illiquid shock. 

Corollary 3.7. By Proposition 13.21 the price under the minimal entropy martingale mea- 
sure Q E , denoted p E (t,S) for the case M t = 0, is given by substituting i>oi with Dqi, and 
l'io with z>io in formula ()3,2|) . 



4. Exponential Indifference Pricing 

To account for investor risk preferences, we next investigate a nonlinear pricing mecha- 
nism in the presence of liquidity shocks. In this section, we derive the exponential indiffer- 
ence price of /i(St) using methods of stochastic control and nonlinear partial differential 
equations. Throughout this analysis, the exponential utility function is used: 

u(x) = -e" 7X , (4.1) 

where 7 > is the investor's risk aversion parameter. Exponential utility has various 
advantages, such as yielding prices that are wealth independent, and being related via the 
dual representation to the minimal entropy martingale measure. 

A starting point for the approach of utility-based pricing are the value functions U l (t, X, S) 
for the optimal investment problem in the presence of a contingent claim, 

U l (t,X,S) = sup El X)S:i \- e -^ XT+h ^] , i = 0,1. (4.2) 
Standard stochastic control methods imply the following: 

Proposition 4.1. The optimal value functions U^t, X, S) with (t, X, S) G [0, T] x R x R+ 
for a holder of a bounded contingent claim paying off H{St) at terminal time T are given 
by 

U\t,X,S) = _ e -7* e -7^,S) ; 

where R l (t,S) are the unique viscosity solutions of the coupled semi-linear system 

R° t + -a 2 S 2 R° ss - ^VtC* 1 -* ) + (jg + ^ = , 

2 T T (4.3) 

7 7 
with the terminal conditions R l (T, S) = h(S), i = 0, 1. 

As a Corollary, taking h(S) = 0, we recover the value functions V l (t,X), i = 0,1, for 
the Merton problem of optimal investment in the liquidity switching market. 

Corollary 4.2. The value function for optimal investment subject to illiquidity shocks is 

V*(t,X)± sup E^[-e^] = -e-^W, i = 0,l, (4.4) 

where Fo(t) and Fi(t) are given in ([3 
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Proof. This is a special case of Proposition 14.11 with zero options. We observe that if 
the terminal condition is zero, the solutions of (|4.3p are independent of (St), R l = K l (t). 
Substituting Fi(t) = e _TRI W and simplifying then yields (|3.4p which are easily seen to be 
smooth. The link between the Merton problem and relative entropy minimization that was 
originally used to derive i^'s in (|3.4p is well-known, see e.g. |15| Theorem 11]. □ 

Remark 4.3. Recall that without liquidity shocks, i.e., z^oi = 0, the Merton value function is 
V(t,X) = -e~ lX F 2 (t). The possibility of illiquidity increases the risk associated with the 
investment behavior. As one would expect, the discount rate has to be therefore increased 
relative to the perfectly liquid case, > F 2 , see Remark 13.51 

4.1. Iterative Approximation. Since we are unable to provide an a priori regularity 
estimate on R l in ()4.3|) . we instead pursue a decoupling approach. The main idea is to 
approximate R° and R 1 through a fixed-point iterative scheme yielding smooth approxi- 
mations i = 0,1, k = 1,2,.... We show that R&Q converge to R l uniformly on 
compacts, and provide a characterization of R^' k ^ as a unique classical solution of the cor- 
responding (de-coupled) HJB equations. The following Theorem summarizes our results, 
with the proof in IB] 

Theorem 4.4. There exist functions U^' k '(t, X, S), i = 0, 1, k = 0, 1, . . . such that: 

(1) (t, X,S) € C^ 2 ' 2 ([0,T] xKxl + ) are classical solutions of the HJB equations 
dp}; 

(2) U^ i,k > (t, X, S) admit the control representation 

m k \t,X,S)= sup El XAi \u(X ( i' k) +h(S$> k) ))\, i = 0,l; (4.5) 

(3) For any compact subsets K x C K and K$ C (0, oo) we have 

lim sup \U ( - i ' k \t,X,S)-U i (t,X,S)\=0. (4.6) 

fe ^°° (t,X,S)e[0,T\xK x xK s 

4.2. Indifference Prices. Using U l and V 1 , the buyer's indifference prices p (initial state 
0) and q (initial state 1) are defined via 

U°(t,X-p,S) = V°(t,X), (4.7) 

U 1 (t,X-q,S) = V 1 (t,X). 

As is well known, under exponential utility the indifference price is wealth independent, 
p = p(t, S), and q = q(t, S). We know the value functions V % = — e~ lX F{(t) from (|4.4p and 
?7*(t, .X", S) = — e lX e~^ Rl<yt,s ^ from Theorem 14.41 so combining these with the indifference 
price definition and simplifying equation (14.3p . we obtain the following PDEs for p and q: 



Proposition 4.5. The option buyer's prices are the unique viscosity solutions of the semi- 
linear elliptic PDEs 

. 1 2 C 2 ^01-^1 -t(o-p) , (^0 + ^0l) l-^o n 

7 t\ 7 7 b\ 

with terminal conditions p(T, S) = q(T, S) = h(S). 

We note that all the terms involving Fi(t) in ()4.8j) are bounded on [0, T\. 



(4.8) 
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Remark 4.6. One can similarly consider the writer's indifference price p w (t, S) and q w (t, S), 
which then solve 



(4.9) 



pf + l -alS 2 p w ss + M^k e 7tt»-P-) _ (*±f*i) + in = 0) 
2 7 F 7 7F0 

W + ^10 ^O e7 (p»-g») _ ^10 + 1^1 = Q 

7 f 1 7 7 i*\ 

with terminal condition: p w (T,S) = q w (T,S) = h(S). 

The utility indifference pricing procedure also provides an explicit identification of the 
hedging position, which is found naturally as part of the optimization problem, namely 
through a maximizer n*(t, S) of (|4.2j) . In an incomplete market perfect hedging is not pos- 
sible, but this utility indifference hedging strategy provides a dynamic optimal adjustment 
of the portfolio strategy for an investor who holds the contingent claim. For exponential 
utility, the optimal trading strategy n*(t, S) satisfies 



Mo _ q 9p(t, S ) 



n*(t,S)X = ^-S^>, (4.10) 



in which the first term is exactly from the classical Merton problem, and the second term 
dP Qg S ^ is the counterpart of the classical Delta to account for hedging motives of h(Sx) by 
the risk-averse investor. We discuss (|4.10p in Section [7.31 

Remark 4.7 (Duality Representation). Recall relative entropy defined in (13. 5|) . One obtains 
(see equation (4.6) in [3]) the following indifference price expression: 

p(t, S; 7) = inf {e^ [h(S T )} + - (H(Q\¥) - H(Q E \¥)) ) . (4.11) 

The representation (|4.1ip shows that the indifference price is decreasing in 7, and yields 
the risk aversion asymptotics (see Proposition (1.3.4) in [T]) : 

lim p(t, S; 7) = inf E Q \h(S T )] , (4.12) 

lim pit, S; 7) = E QB \h(S T )} ■ (4.13) 
7^0 

This implies that when the risk aversion vanishes, the indifference price converges to the 
MEMM price. When the risk aversion goes to infinity (total risk aversion), we obtain the 
super-replication price. 

Remark 4.8. Continuing Remark 12.11 one could consider indifference pricing of /i(S'r) as- 
suming more general dynamics in the illiquid regime. However as mentioned before, this 
would require adding 77 as a state variable when M t = 1 (i.e. for the g-value function in 
(14. 8p ). We present details of the corresponding generalizations of (|4.8p inlDl 



5. Asymptotic Analysis 

The indifference prices from (|4.8p do not have a closed-form analytical solution. To gain 
further insight, in this section we perform regular perturbations of (|4.8p in terms of a small 
risk aversion parameter 7 or a small probability of liquidity shock uq\ . 
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5.1. Small Risk Aversion Asymptotics. For the asymptotic case 7 — >• 0, we can reduce 
the non-linear reaction-diffusion indifference price equation (|4.8p to a linear one, obtaining 
a closed-form representation for the limiting price. 
We apply the formal asymptotic expansion 

p(t, S) = p(°) (t, S) + 7 P (1) (t,S) + ... and q(t, S) = q^ (t, 5) + 7 </ (1) (t, 5) + . . . . 

Theorem 5.1. Lei Bij(s,t) = exp(— Uij(u) du) fori = 0, 1 and j = The expansion 

terms p^°\ p^ 1 ', q^°\ and admit the stochastic representations 

fT r 1 

dr, 

(5.1) 



p<°> (t, S) = B 01 {t, T)P BS {T -t,S) + J % (t)B 01 (t, r)E tjS [q^ (r, S T ) 
q^(t, S) = B w (t, T)h(S) + ^ u w (r)B w (t, r)p<°) (r, S) dr, 



\( q {0) (T,S T )- P ^(T,S T )) 2 -qV(T,S T ) 



P (1) (t,S) = -J U 01 (T)B 01 (t,T)E t 

q {1) (t,S) = -j\ (r)B 10 (t,r) Q {p^ (r, S) - q^ (r, S)} -pW(r,S) ) 



dr, 



(5.2) 

where Pqi (i) P10 (i) are specified in Proposition \3.b\ see (|3.9j) . 



Proof. Formally matching terms in powers of 7 in the PDE (|4.8p yields the following 
equations: 

P f (0) + ^0 2 ^S + %(i)(<? (0) - P (0) ) = 0, (5.3a) 

+ M*)(P (O) -? (O) ) = 0, (5.3b) 

Pf } + \^ P f s + v Ql {t){qU -p«) - ^ 01 (*)(s (0) -pW) 2 = 0, (5.3c) 

g« + ^ (t)(p (1) " - ^io(*)(<7 (0) - P {0) ? = 0, (5.3d) 

with boundary conditions: 

p(°)(r,5) = g (0) (T,S) =fc(S), and p (1) (T, 5) = g (1) (T, 5) = 0. 

The equations for p(°) and p*- 1 ) are typical Cauchy problems described in \i'6\ Ch. 5.7]. 
Since we have constant diffusion and volatility terms, and the source and the terminal con- 
dition functions are bounded, there exists a unique solution to (|5.3p by a direct application 
of the Feynman-Kac representation in Theorem 7.6 in [13j . leading to (|5.ip - (|5.2|) . Finally, 
the equations for q(°> and q^ are linear first-order ODEs that have an explicit solution 
using an integrating factor. □ 

Using the dual representation (|4.12p . it is immediate that p(°) = p E and g(°) are the 
MEMM prices, see Corollary [577] and pp. 

5.2. Small Probability of Liquidity Shocks. Under normal conditions, the likelihood 
of a liquidity shock is small (typically caused by a natural or geopolitical catastrophe that 
is rare by definition). In terms of our model, the transition probability i^oi is small under 
this scenario and we may view the resulting option price as a perturbation of the classical 
Black-Scholes price corresponding to uq\ = 0. In particular, for a given option maturity T, 
the probability of more than one liquidity shock is exceedingly small, O(^qi)- 
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In light of these considerations, we analyze the asymptotic case of at most a single 
liquidity shock. Such a restricted model may be obtained by redefining the state-space of 
{Mi ) as E = {0, 1,2}. As before, we assume the market starts with the regular liquid state 
(0). Due to the unpredictable illiquidity shock, the market is susceptible to transiting to the 
illiquid state (1), and before the terminal time T it may come back to the absorbing liquid 
state (2) which is now immunized from any further illiquidity shocks. The corresponding 
infinitesimal generator is 



A 




(5.4) 



Similarly, we define the regime-dependent discounting factors to the utility functions in 
Merton's problem, denoted by Fq, F\ and F2. 

-do(T-t) 



F 2 (t) = F 2 (t) = e 



F 1 (t) = F 2 (t) 
F (t) = F 1 (t)^ 



d 



^10 



d 



^10 

do 



do + i>oi - ^10 
which are the solutions of the ODE system 

F'(t) = (D-A)F(t), F(T) = 



d - vio 

-(d +u 01 )(T-t) _ 



-u 10 (T-t) 



(5.5a) 
(5.5b) 

(5.5c) 



with D 



( 





o\ 


2^ 















1° 





Mo 

2^J 



Note that F 2 (t) = F 2 (t), the discount factor in the classical Merton's problem. This is 
not a surprise considering the state (2) in this three-state liquidity switching case is in fact 
the classical complete market. In addition, we have -F (i) < and F 2 (t) < Fi(t). We 

obtain that the MEMM in the single-shock market above is given by Q E = Q d , where the 
minimal entropy risk factor is a(i,j) = "p^y> hj S {0, 1, 2}. Thus, under Q E the Markov 
chain (Mt) has transition intensities 



^01 (t) 



Let us denote Bij(s,t) = e 



4 p~ fs i> ij( u ) du 



and ^io(^) - 
= 0, 1 and j = 1 



v P2(t) 

V10 WY 



(5.6) 



Corollary 5.2. Using the pricing formula in Proposition 13.21 the MEMM price p E of a 
contingent claim with payoff h{Sx) in the single-shock market is 



p E (t,S) 



L(t) 



B 01 (t,T)h(S T ) + / £>oi(T)B 01 (t,T)L(T)dT 



where 



i>io(s)Si (r, s)h(S T - s + T ) ds + B w (t, T)h{S T ). 



(5.7) 
(5.8) 



Next, the corresponding indifference buyer's price assuming initial regime Mt = is the 
unique classical solution of the semi-linear PDE 



Pt + PSS 



vqi g(t,S) 

7 F 



+ 



(do + m) i(Fi 



OH 



7 ^o 



0, 



(5.9) 
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with terminal condition p(T, S) = h(S) and the explicit source term 

g( t , S) = t Vwe ^u e ~d {T-t-u) e -^P BS {T-t~u,S) du + e - Vl0 (T-t) e -yh(S)_ 



As in last subsection, we can further expand in 7 to obtain p(t,S) = p(°\t,S) + 
7P (1) (t,5) + . . ., with 

(t, S) = B01 (t, T)P BS (T -t,S) + ^ voi (r)Boi (t, r) 

x l^^oCO^ioCr.OPBsCr-C.^dC + BioCr.rjPBsCr-t.^ldr (5.10) 

P {l) (t,S) = ~-Et, s \ F *ai(T)8 01 (t,T) x { f T v w (()B 10 (t,0 {p {0) (r,s T ) 



2 

, 5c [M5r)]) 2 dC + B w (T,T)(pW(T, S T ) - h(S T )) 2 }dT] . (5.11) 

Compared to the original model we not only have a single semi-linear PDE for the 
indifference price p but also explicit expansion terms and p^\ The zero-order term 
p(°) can be expressed as the weighted average of the Black-Scholes prices of options with 
same initial condition but varying maturities, 

p(°\t,S)= [ P BS (T-u,S)w(du) (5.12) 



for a probability measure w(-) on [t,T] defined implicitly in ()5.10j) via the transition in- 
tensities Vij(t). Similarly, p^ is clearly negative and is a weighted average of squared 
differences p(°\u, S) — Pbs{T — u, S) with same weights w(du) on [t, T] as in fl5. 12f) . 



6. Implied Time-to-Maturity 

In our liquidity switching model, the random trading interruption is the only extra risk 
in comparison to the Black-Scholes model. As the market is static during the no-trading 
time periods, one could naturally regard the model as a stochastically shortened time-to- 
maturity Black-Scholes model. In this subsection, we accordingly define the concept of 
implied time-to-maturity as the measure of the investor's view of effective time horizon 
based on her pricing approach. Let us first define the realized time-to-maturity T, 

T= / l{M e =o}ds. (6.1) 
J o 

Then T is a bounded (/^-measurable random variable. Due to the independence of (Wt) 
and (Mi) and using the strong Markov property of Brownian motion, we can view the 
stock as following (|2.3p . and the contingent claim as paying h(Sj-) at the random terminal 
date T ■ 

Definition 6.1. The Adjusted TTM T under a measure Q is defined as 

T*(T; <Q>) = ~E Q [T\M = i], i = 0, 1, (6.2) 

i.e., the original maturity T subtracted by the expected no-trading time duration. 

Assuming a time-stationary generator A as in (I2.ip . standard analysis of the sojourn 
times of the 2-state Markov chain (Mi) yields 



EUROPEAN OPTION PRICING WITH LIQUIDITY SHOCKS 



12 



Corollary 6.2. The adjusted time-to-maturity under the minimal EMM Q° is given ex- 
plicitly by 

T°(T;Q ) = i ^ Ui + vio(m + vio)T - z, ie-^ 1+ ^ T j , 

l n V 6 - 3 ) 



T 1 (T; i 



(^oi + ^io) 2 - 



-^io + v 1q {vqi + viq)T + v w e 



-(v 01 +v 10 )T 



Simple algebra arguments lead to < T l (T) < T°(T) < T. The adjusted TTM simply 
subtracts the expected period of market freeze, without regard for when it will occur. 
Thus, it does not account for option-specific features whereby the timing of the shock 
is also important, as the option time-decay depends on t. This motivates the following 
definition (compare to the popular concept of implied volatility). 

Definition 6.3. Given model-based price p = p(t,S), the implied time-to-maturity 

T(t,S;p) is defined implicitly via 

P BS (f(t,S),S)±p(t,S). (6.4) 

In other words, T is defined as the time-to-maturity that equalizes p with the classical 
Black-Scholes price. 



Remark 6.4. Certainly this implied TTM concept is meaningful only when (|6.4p has a 
unique solution, i.e. when the time value of an option is meaningful, such as for Calls and 
Puts. In such cases, we could use T to parameterize model option prices. However, for 
options, such as digital Calls, whose Black-Scholes price is not monotone in T, there is no 
unique implied TTM and thus this concept should not be used. 

The implied TTM can be used to explain the spread between the Black-Scholes price 
and the liquidity-adjusted valuation. Recall that the option value may be decomposed 
into intrinsic value (payment an agent would gain if the option is exercised immediately) 
and time value (the value of delaying exercise until maturity) . Typically, the time value is 
increasing in time-to-maturity. Consequently, since a potential liquidity shock shrinks the 
expected TTM, the option price decreases (note that the intrinsic value is unaffected due 
to our assumption that (St) is static during the shock) due to the extra time-decay. 



7. Numerical Examples 

In this section we present a series of numerical experiments to illustrate our derivations 
in the liquidity switching model. Recall that we have obtained three main types of prices, 
the indifference price p from (j4.8|) , the MEMM price p E in Theorem 13.71 (equivalent to p^ 
in Theorem 15. ip , and the minimal martingale measure price p MM . The indifference price 
is obtained through utility maximization arguments, while the other two are from special 
martingale measures. Both p E and p MM are easily computed from the generic formula 
(|3.ip using the explicit descriptions of the corresponding transition intensities of (Mt). We 
compute p by solving the corresponding semi-linear PDE using a mixed implicit-explicit 
finite difference scheme, see ICl for details. 

Table [U summarizes the parameter values used. Note that with z/qi = 1 and z^io = 12 we 
assume that liquidity shocks occur at a rate of once per year and last an average of one 
month. As two representative types of European options we consider vanilla Calls with 
payoff }i(St) = (St — K)+ and digital Calls with payoff h(Sx) = l{s T >K}- 
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Parameter 


Meaning 


Value 


Mo 


Growth rate of the stock 


0.06 




Volatility of the stock 


0.3 


^01 


Transition rate from state - 


+ 1 1 


^10 


Transition rate from state 1 - 


+ 12 


K 


Strike Price 


10 


T 


Option's maturity in years 


1 


7 


Risk aversion parameter 


1 



7.1. Price Comparison. We first discuss the observed relations between the different 
model prices p, p E , and p MM . As the MMM price p MM assigns zero illiquidity risk pre- 
mium, we have p E < p MM } see Proposition 13.61 Also due to risk aversion, p < p E , so that 
we have 

p < p E < p MM . 

Numerically, we find that at the given frequency of liquidity shocks, the difference between 
the MMM and MEMM prices, p MM vs. p E , is negligible across all option types we tried, 
showing that the entropic liquidity risk premium is generally very small. Next, our nu- 
merical results confirm that p E > p. However, the magnitude of the price spread varies 
for different types of options. For Call options, there is no significant difference (less than 
0.5% in all cases we tried); while significant spread is found for digital options, see Figure 
[TJ In absolute terms, the spreads are largest at-the- money (right panel of Figure [TJ, how- 
ever since the Call price is increasing with respect to S, the percent difference is generally 
decreasing in S (left panel of Figure [TJ . 

Regarding the price spreads versus time-to-maturity (TTM), intuitively, as TTM shrinks, 
there is less chance to have liquidity shocks, thus smaller spreads should be expected. This 
pattern is observed for Call options. However, for digital Calls, the effect of TTM on p E —p 
is ambiguous. This is due to the digital Call's large Theta (JjO around the strike. Hence, 
for short-maturity ATM digital Calls a large liquidity cost is observed, see the hump at 
t = 6/12 in the right panel of Figured) Moreover, out-of-the-money, the digital Call B-S 
price is not monotone in time-to-maturity. Therefore, the shorter TTM (t = 6/12), the 
smaller p E , and thus the highest percent spread happens at the deep OTM case, see left 
panel of Figure [TJ 

Remark 7.1. In this paper we have chosen a probabilistic description of liquidity shocks 
through the Markov chain (Mt). An alternative formulation could focus on worst-case 
analysis, treating the shock interval [toi,tio] as unknown constants ,see [18]. In such min- 
max setting, the liquidity spread could be large even as t — > T if the option time-decay is 
strong. 

Next, we study the quality of the asymptotic expansion in vqi in Section [5T2"1 bv comparing 
the indifference option price p(t,S) and its approximation p(t,S) from (15. 9|) . Recall that 
p(t, S) corresponds to the price assuming at most one liquidity shock on [0, T] and is 
therefore the first-order term in the asymptotics using vq\. We find that this first-order 
approximation explains over 50% of the spread across most stock prices (see Figure [2j and 
can therefore provide a useful rule-of-thumb without the need to solve any nonlinear PDEs. 

To understand the nonlinear nature of the indifference pricing rule, Table [2] summarizes 
the per-contract prices as we vary the quantity of contracts to purchase/sell. As expected, 
the indifference price decreases as the number of contract grows due to the increasing risk 
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0.15 



0.05 









f = 6/12 











10 12 14 

Stock Price S 




10 12 14 

Stock Price S 



Figure 1. Price difference between MEMM p E (t, S) and indifference price 
p(t, S) in percent (left panel) and in dollars (right panel) for buying 10 
digital Calls. Parameters are from Table [TJ 




10 12 

Stock Price S 



Figure 2. Small foi asymptotics. We plot the percent price difference 
between p(t, S) defined in (|5.9p and p(t, S) for buying 10 digital Calls. All 
the option and model parameters are from Table [TJ 



of holding more options. Also, the linear risk neutral prices (MEMM, MMM and adjusted 
B-S) are between the buyer and writer's prices for n = ±1 (recall that MEMM is the 
limiting per-contract price as n — >• 0). Recall that with exponential utility contract volume 
and risk-aversion asymptotics are equivalent, so Table [2] can also be interpreted in terms 
of modifying the risk-aversion of the investor. 

7.2. Implied Time To Maturity. As we discussed in Section [61 it is tempting to charac- 
terize the model prices using the Black-Scholes equivalent (i.e. implied) time-to-maturity 
defined in (16. 4p . As the Black-Scholes digital option price is not a monotone function of 
TTM, it will only make sense to discuss Call/Put options' implied TTM. The left panel of 
Figure [3] shows the adjusted T°(T - 1; <Q°) and implied TTM f(t, S) as we vary the overall 
time-to-maturity T — t. We observe that T°(T — t) is almost linear in t and provides a 
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Table 2. Prices per option contract. Positive number of contracts n refers 
to buyer's prices and negative n to writer's prices. The parameters used are 
from Table [TJ 



Indifference Price per Contract 



# of Contracts n 


Call Option 


Digital Call Option 


in Possession 


So — 8 


S = 10 


So = 12 


5*0 — 8 


Sb = 10 


S = 12 


10 


0.2875 


1.0720 


2.4476 


0.1655 


0.4229 


0.6705 


5 


0.3128 


1.1264 


2.4872 


0.1751 


0.4370 


0.6826 


1 


0.3222 


1.1442 


2.5014 


0.1793 


0.4433 


0.6883 


-1 


0.3253 


1.1496 


2.5060 


0.1811 


0.4461 


0.6909 


-5 


0.3296 


1.1573 


2.5124 


0.1856 


0.4535 


0.6980 


-10 


0.3333 


1.1635 


2.5178 


0.1967 


0.4723 


0.7155 


Risk Neutral Price 


MEMM p*(0,S) 


0.3235 


1.1466 


2.5034 


0.1801 


0.4447 


0.6897 


MMM p MM {0,S) 


0.3236 


1.1467 


2.5035 


0.1801 


0.4447 


0.6897 


AdjB-S P BS (T°(X),S) 


0.3247 


1.1496 


2.5053 


0.1798 


0.4425 


0.6865 


B-S P B .s(T,S) 


0.3534 


1.1924 


2.5441 


0.1857 


0.4404 


0.6764 



reasonable approximation to the indifference TTM T(t, •). This implies that using an ad- 
justed TTM within the classical Black-Scholes formula can be used as a simple correction 
to take into account liquidity risk. 

The adjusted TTM T is generally larger than the implied T due to the risk aversion and 
non-constant option time-decay. However, we find that this relationship is not universal; 
in fact for deep ITM or OTM options it may be reversed. This is illustrated in the right 
panel of Figure [3] that shows dependence of the implied TTM on the initial stock price So- 
We observe a "smile" shape such that the most time- value lost is for at-the-money options. 
This is consistent with intuition that liquidity risk is related to the Charm and Gamma 
greeks which are largest at-the-money. 



■(T-t)-T(t,S) 
-(T -t) - f°(T -t) 




10 12 14 16 

Stuck Price 5 



Figure 3. Adjusted/implied TTMs for a European Call. The parameters 
are as shown in Table [TJ In the left panel we take So = 10; in the right-panel 
we have t = 0. 



7.3. Hedging. Recall that in the problem defining U l (t, X, S), the optimal stock holding in 
dollars is ir* = — S dp g^ . The first term is the constant arising in the Merton's hedging 
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problem, and the second term is the Delta hedge analogous to the classical Black-Scholes 
case except the price is the indifference price p(t,S). Figure H] compares this indifference 
hedge to the classical B-S Delta Abs(T, S) = l9Pg gg T '' 5 '^ using the original time-to-maturity 
T, or the adjusted time-to-maturity T°(T — t; Q°) from (|6.3D . 




Figure 4. The Deltas A BS {T - t,S), A B5 (T°(T - t),S) (adjusted B-S 
delta) and dP Qg S ^ (indifference price delta) at t = for buying 10 digital 
Calls with strike K = 10. We compute the sensitivity ^ directly from the 
PDE (14. 8p using a first-order finite-difference approximation. The parame- 
ters are as shown in Table [H 

Using Abs to denote the classical Delta in the Black-Scholes model, and recalling that 
p(t,S) =: PBs(T(t, S), S) we decompose 

d P (t, S) 



dS 



A BS (T -t,S) + A B s(T"(T -t),S)- A BS (T - t, S) 



adjusted TTM spread 

+ A BS (f (t, S), S) - A BS (f°(T - t), S) + - A BS (f(t, S),S). 

implied TTM spread TTM smile correction 

The three terms in the above comparison between the classical Black-Scholes Delta and the 
indifference hedge correspond to adjusting for expected time- value loss, further correcting 
for timing of the shocks and the risk-aversion (by using implied TTM T) and finally a 
remainder term due to the fact that T(t,S) itself depends on S. Using Abs(T\, S) — 
Abs(T2, S) ~ d^-g^i'^) (j^ — T2), we may then relate the Delta spread to the joint effect 
of the charm (or Delta decay, measuring the instantaneous rate of change of delta over 
the passage of time) and the TTM spread. Since we have that T°(T — t) < T — t and the 
digital Call charm is positive unless the option is deep ITM/deep OTM, we have 

A BS (f, S) > A BS (T, S), S/K e [0.7, 1.3]. 

The sign of dP gg S ^ — Abs(T, S) is more complicated, since implied TTM is ill-defined for 
digital options due to the changing sign of Theta. Overall, we observe that liquidity shocks 
make the hedges more extreme; thus for OTM options, the risk-averse investor tends to 
hold fewer shares of stock; for ITM options the investor tends to hold more. An intuitive 
understanding is that the illiquidity shock reduces variance of St, so starting in-the-money 
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one is more likely to end ITM, and hence the investor should hold more shares there, and 
vice versa. 

8. Conclusion 

In this paper, we considered the valuation and hedging problem of European options in 
a regime-switching market subject to liquidity shocks. Investors exposed to such illiquid- 
ity risk will generally demand a positive risk premium (i.e. offer to pay less to purchase 
options). Our numerical experiments suggest that the entropic risk premium associated 
with the MEMM Q E is negligible in our model. On the other hand, indifference prices may 
impose significant liquidity premia depending on the option type. While these are small for 
vanilla Calls and Puts, we observed up to 10% price differences for digital options (namely 
OTM digital Call). For options with monotone time- value a simple way to adjust for liq- 
uidity risk is available through the adjusted TTM formalism which allows continued use of 
the Black-Scholes formulas after modifying the option's time-to-maturity. We believe this 
idea is a useful guide for traders who want a transparent rule to adjust their valuations 
from the base-case. In terms of the resulting hedging adjustments, the illiquidity threat 
causes the investor to hold more shares in-the-money and fewer shares out-of-the-money 
due to the reduced time-to-maturity; this correction can be approximated using adjusted 
Delta as well, and more generally by looking at the resulting TTM impact and the option 
Charm. 

Even after several simplifying assumptions, our model remains complex. In|D]we revisit 
a more general formulation of the impact of liquidity on asset dynamics; as can be seen in 
such generality little can be said beyond numerical studies. Consequently, we recommend 
to start with a well-understood model (e.g. Black-Scholes) and then incrementally adjust 
for other risks (such as risk aversion, jump risk, illiquidity risk, etc.); in that case our 
analysis here, isolating impact of trading interruptions, can be used as one of such steps in 
obtaining the final risk-adjusted option value. 

Appendix A. Proof of Proposition 13.21 

Under a general EMM Q Q , at time t, we prove that the risk neutral price p a (t, S, 0) 
with initial state St = S, Mt = can be written in the form (13. 2|) . Similar arguments for 
p a (t, S, 1) will be omitted. 

Let Vij = at(i,j)uij be the transition rates of (M t ) under Q a . Denote by roi (resp. 7~io) 
the first transitioning time from liquid state to illiquid state 1 (resp. from state 1 to 
state 0), which have intensities uqi and viq. Therefore, starting with Mt = 0, we have 
two possible cases, represented by {tqi > T}, and {t < tq\ < T} leading to the following 
expansion of the price 



1 {t<T i<T}\ 

= Q^ i0 (roi > T)E« S [h(S T )) + jT goi(T)E« STtl [h(S T )] dr 

where the last equality is due to the independence of (St) and (iVoi(t)) and goi(-) is the 
density of Tq±. Substituting Q a (roi > T) = exp(— J*. T i>Qi(u)du), and 

d 



9oi (t) = ^Q?(r i < s) ^ = %l(t) exp(- J D 01 (u)duj , 
yields (1331) . 
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Appendix B. Proof of Proposition 14.11 and Theorem 14.41 

Let us first sketch out the main four steps of the proof. First, through the standard 
stochastic control arguments, we derive the following HJB equation for U l (t,X, S), 

vf + fi su° s + 1 -<j 2 s 2 u° ss + miu 1 - U°) (B.l) 

+ sup^o^Tr + \olX 2 U Q x y + a 2 SXU° Sx n} = 0, 
U} + v w (U° -U l ) = 0, 

with terminal condition U l (T, X, S) = — e~~ / ( x+h ( s ^ . Using the scaling properties of expo- 
nential utility, we consider a candidate solution of the form U l (t, X, S) = —e lX e~' yRl<yt,s \ 
Substituting into (jB.ip yields that R l (t,S) must satisfy the equation (|4.3|) . The difficulty 
arises due to the coupled form of (|4.3p which poses several technical challenges. 

To circumvent this problem, in our second step, we consider an approximation scheme 
i?(*' fc ) defined in (|B.2p for W and show that the corresponding functions R^' k ) are bounded 
and are classical solutions of the respective PDEs. This is the longest part of the proof and 
relies on a doubly-iterative setup. Third, we provide a stochastic control representation 
for which shows the convergence Hindoo -> R\ Finally, using the properties 

of we establish a verification theorem for R l . 

To introduce our iterative scheme, we consider a model with at most k £ N liquidity 
shocks on [0,T] (a similar strategy is used in j5]). Formally, let 

K{t) := max{k : Tqi < t} 

denote the number of liquidity shocks by time t. Then we define a counting process (N^t)) 
with intensity ^oil{Af t =o}I{ftT(t)<fc} an d use the superscript k to denote the resulting wealth 
process (X-) and stock process (S k ). In analogue to the original U l (t, X, S) we then 
consider the value functions U^ l ' k \t,X,S) defined in (|4.5p . 

Step 1: Control Representation. Usual arguments imply that U^ l ' k ) from fj4.5|) satisfy 
(in viscosity sense) 

£/(0,0) = _ e -d (T-t) e - 1 X e - 1 t[h(S T )}^ 

Of*) + u w (U^ - U&Q) = 0, k = 0, 1, . . . 

(j<U<) + lotSPugf + v SUp k) + voiiU^-V - uiW) 

+ smv{\<jIX 2 U^k 2 + + hqXSU^tt} =0, k = 1, . . . , (B.2) 

where the terminal condition is U^' k \T,X, S) = - e ~^ x e~^ h{ - s \ i = 0,1. 
Making the ansatz 

U^ k \t,X,S) = -e-^ x e-^ R( °- k) ^ and U^ k \t, X, S) = -e^ x R^ k \t, S) (B.3) 
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and substituting the candidate optimizer Tr^'°^'*(t, X, S) = —-^r^— 2 - — . ( ,k) XS — into equation 
(|B.2p . we are led to consider the following iterative equations for R^ l ' k ^s, 

(B.4a) 

R^k)(T,S) = e-^ S \ (B.4b) 
0, R {0 ' k) {T,S) = h{S), (B.4c) 

where k = 0, 1, . . . and 

S , Q, R) 4 ^Qe-r* - *±M. (B.5) 
7 7 

Once we show in Step 2 that there exists a classical solution to (|B,4cp . our ansatz (|B.3P 
will be validated and the control representation (|4,5p immediately follows. 

Step 2: Classical Solutions. The present step is to prove that equations (1B.4|) have 
bounded classical smooth solutions R^ l ' k \t, SYs for t € [0, T], € (0, oo), i = 0, 1, and 
< oo. 

By assumption, h(S) is bounded, so there exist constants A, B such that A < h(S) < B, 
for all S. It follows that A < P B s(T -t,S) = E[h(S T )\ < B, for all t and S. Being related 
to the classical B-S price, R^°\t, S) is in C^ 2 ([0,T) x R + ) and satisfies 

A < R(°'°\t,S) < B + — (T — i). 

7 

Next, i?( 1,0 ) follows a linear ODE in (IB,4bj) . thus we have explicitly 

R(W)(t, 5) = e -7MS) e -io(T-t) + ^ T ^ oe -7^ '(«.s) e -^(«-*) dti . 

It is easy to see that e -yB-MT-t) < fl(i.o)( 4j 5) < e -7^ y(t, 5). 

We now use induction to show that there exists a series of smooth and bounded functions 
RP'Q and R&Q, k = 1,2,..., which are the solutions to equations (jB.4|) . To apply the 
inductive argument it is sufficient to show that there exists a classical solution R^ ' 1 ^ to 
the semi-linear parabolic PDE ([Hie]) such that A < R^(t,S) < B + &(T - t) (same 

argument as for it^ 1,0 ) will then take care of R^ 1 ' 1 ^). Typical existence and uniqueness 
results for solutions to semi-linear parabolic equations, such as [T], do not apply here 
because they require that the non-linear term f(t, S, Q, R) be uniformly Holder-continuous 
on the domain. In our equation, however, the non-linear function f(t,S,Q,R) is not 
globally Holder continuous in variable R. We therefore follow Zhou [26] to resolve this 
difficulty through a further iterative argument. The main idea is to convert (|B.4cP into 
a series of linear PDEs, and using the comparison principle prove that the linear PDEs 
solutions converge to the solution of the original semi-linear PDE. 

First, we prove that if the solution R^ 0,1 '(t : S) exists, then it is bounded. Let us define 

R(t,S) = B + — (T -t) and R(t,S)=A. 



R^°Ht, S) = P BS (T - t ,S) + -±(T-t), 

7 

R^-u 10 R^ + u 10 e-^ k) =0, 

R (W + la 2 S 2 R { ^ k) - /(i.^i?^- 1 ),^)) = 
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Then 

R t + ±<j 2 S 2 R ss = -^ <f(t,S,R^°\R), R(T,S) = B>h(S), (B.6) 
2 7 

Rt + l^ S 2 R ss = 0>f(t,S,R^°\R), R(T,S)=A<h(S). (B.7) 

So by the comparison principle stated in section (i) in the Appendix of |26j . we have 
A < Rp'Vtt, S) < B + &(T - t) if it exists. Note that f(t, S, Q, •) is not globally Holder 
continuous, but it is locally Lipschitz continuous in R when R £ [A,B + -^T]. Therefore, 
there exist positive constants C_ and C such that 

-C(R} - R 2 ) < f(t, S, Q, R l ) - f(t, S, Q, R 2 ) < C(R l - R 2 ) (B.8) 

for A = R < R 2 < R 1 < R < B + ^-T and S G M + . Note that C and C are independent 
of R 1 , R 2 , S and t. 

Next, we construct the sequence JRC'Ww) via #( >i)(°) = R and for j = 1, 2, . . ., 

R (01M + CR mU) = fa S , RWV-V), where (B.9) 
CR= -a 2 S 2 R S s-CR, and f{t,S,R) = f(t,S,R^°\R)-CR, 

with terminal condition R^'^^^T, S) = h(S). The PDE (|B.9|) is linear, so classical exis- 
tence and uniqueness results guarantee we have a series of smooth solutions R^ 0,1 ^',j = 
1,2,.... Next we use induction to prove that R < ••• < #(°>i)tj) < #(o,i)0'-i) < ... < 
fl(0,l)(l) < fl(0,l)(0) = #. 

First, we prove i? < fl^X 1 ) < J R(°. 1 )(°) = #. From (lR6l) and l|R9|) . we have 

(#-#(0:l)(i)) f + £(#_#(o,i)(i)) <0, 
(fl-i^WXT.S) > 0. 
On the other hand, similarly, we have 
(#(0.1X1) _ R) t + £(#(°.i)(i) -R) = f(t,S,R)-CR + CR 

< f(t,S,R) - f(t,S,R) -CR + CR (by dEB) 

< (by {EE}), 
(i?(°.i)(i))(T,5)- ; R>0. 

Therefore, by the comparison principle, we have R < J R(°'i)( 1 ) < #(°>i)(°) = R, Now we 
assume £ < flC - 1 )^) < #(°>i)0-i) < . . . < # ; and prove that # < R(°^ +1 ) < Rptiti). 
By (|B~9j) . we obtain 

(#(0T)0) _ jj(o.i)(»+i)) t + cfBf- ' 1 ^ - J?^ 1 )^ 1 )) 
= fa,S,RQ> 1 W- 1 ))-f(t,S,RP' 1 M) 

= f(t,S,R^°\R^-V)-CR^-V - \f(t,S,R& \Rt°> 1 M)-VR<- > 1 M] 
= [f(t, S, ^(1,0)^(0,1)0-1)) _ /(tj Sj i2 (i,o) j #(0,1)0))] _ c[#(o,i)0-i) _ #(0,1)0)] 

< (by {EH])), and (i? (0 ' 1)(j) )(T, S) - (J?( ' 1 X fc+1 >)(T ) S') = 0. 
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On the other hand, we have 

(jR (0,D(i+D _ R)t + jC (ij(0,l)0-+l) -R) 

= f(t,S,R^°\RW^) -CRWW + CR 

< [f(t,S,R {1 '°\R^ 1)(j) ) - f(t,S,R (1 ' 0) ,R)} -C[RP> 1 M -R] (by ([BTfj) ) 

< (by (EH])), and (i^N 1 ))^, 5) - R(T, S) > 0. 

So, by the comparison principle again, we have R < j R(°' 1 )0'+ 1 ) < itfMU), completing the 
induction on j. 

Lastly, we prove that the sequence 

{#(0,1) (j)} 

converges to the solution of (|B.4cp for 
k = 1 as j — > oo. Notice that {Rp'W)} is a monotone decreasing sequence, bounded 
below by R, for all (t,S) € [0, T] x M + , hence, it converges to some A 0,1 \t,S) pointwise. 
Using the fact that it^ ' 1 )^) are bounded, it follows that f(t, S, R) is of at most exponential 
growth in S so by Appendix A in [26] we obtain that r^ 0,1 ^ is in C ll2 ([0,T] x M+) and hence 
is the unique classical solution of (|B.4cP for k = 1. This completes the inductive step in k. 

Step 3: Convergence. To establish the uniform convergence in ()4.6p . we apply Theorem 
3.2 in [5] since the exponential utility u(x) = — e _7X is polynomially bounded at 0. It 
follows that for k large enough, U^' k ^ is arbitrarily close to U' 1 which is intuitive since 
the probability of more than k shocks is exponentially small in k. However, since we 
do not have control on the rate of convergence, we cannot establish such properties for 
the corresponding controls ir^°' k ^'* so it is not clear whether the hedging strategies also 
converge. 

Step 4: Verification Theorem. We are finally in position to provide a verification 
theorem for the original PDE (|4.3p . Similar arguments to those below also provide a 
verification theorem for the iterative (|B.2p and are omitted. 

Lemma B.l. Suppose that we have J l (t,S,X) and tp(t,X, S), such that 

(1) f l (t,X,S),i = 0,1 is in C^ 2 ' 2 ([0,T) xlxl+) solve the HJB equations gjfy ; 

(2) JgcroS + J x aoirX £ L 2 for all admissible control laws (irt); 

(3) tjj(t, X, S) is an admissible control; 

(4) For each fixed (t,X,S), 

™p{±a 2 X 2 + fioXJlir + fi XSJ° Sx ir} = ^a 2 X 2 J° xx ^ 2 + (f, XJ° s + f, X S J° Sx )^ . 
Then the optimal value function of (|4.5p is U l (t, X, S) = J l (t,X,S). 

Proof. Given an arbitrary admissible strategy (irt) £ A, we define the process as 
the solution to the equation dX^> = l^ Mt=0 y[fj,()irX^ dt + o^itX^ dWt]- Applying Ito's 
formula to ()B.2|) . we obtain 

dJ°(t, X t W , S t ) = J° dt + JKfioSt dt + a S t dW t ) + \4s4St dt 

+ J%(fi nX t dt + o-QirXt dW t ) + - J^ x a 2 ir 2 X 2 dt + (J 1 - J°) dN 01 (t). (B.10) 

Note: for notational convenience, we omitted the parameters (t,X^\St) in J°. 

Then, in the above equation we replace cL/Voi with d/Voi + voidt, where Nq\ is the 
compensated Poisson process associated with Nq\. Thanks to the integrability assumptions, 
the stochastic differentials in (jB.lOh are stochastic differentials of P-martingales. Therefore, 
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by using the boundary conditions of the HJB equations (|B.2p . integrating up to the first 
liquidity shock instant r^, and taking expectations on both sides under P measure, we 
have 



E: 



t.X,S 



I {T i i<T} J^X^S^) + 1 {tL > t] u{xP + h(S T )) 

i-tL AT 



= At, x, s) + eJ XiS [ jT 01 { j? + \ols> 4 S + h>sj° 8 + l r lx*j» x y 

+ fi XJ s 7r + fi XSJ Sx + ^{J 1 - J°)n}dT . (B.f 1) 

The parameters in J-functions in the integral are (r, X? , S T ) in the right hand side of 
equations. Since the J's satisfy the HJB equations (|B.2p . the term inside the expectation 
on the RHS is negative a.s., and therefore 

Kx,S [hr^JHr^X^S^) + I {T i> r} «(lW + h(S T ))] <J°(t,X,S). (B.12) 

Using the fact that (Xt, St) are static on {Mt = l}i s ° that 

j\t,X,S) =El X)S [l {Th]<T} J (rl ,X,S) + l {T ^ T}U (X + h(S)) , 

and applying the strong Markov property of (Nij(t)) at we obtain 



E 



t,x,s 



l {T ^J°(ri ,X^\S T ^) + I { ^> T} u(X^ + h(Sr)) <J°(t,X,S) 



By induction on the transition times (tqi,^) it follows that 



E 



t,x,s 



Vf < 



T} J°(t? , xg, S t , o ) + l Ko > T} u(X^ + /i(S T ))j < J u (t, X, S). (B.13) 

Since rf — >• oo a.s., and using the fact that J°(-,X. , S.) is bounded from below, we 
conclude using the dominated convergence theorem that 



J u (t,X,S) >El x<s -e W l--f{xP + h(S T )} 



(B.U) 



On the other hand, if we choose the control law ip, the term inside the expectation in 
(jB.lip is zero a.s., yielding 



E- 



t,x,s 



exp(- 7 {4 ,/ ' ) + HSt)})] = J°(t,X,S). 



Combining (jB~T4|) and (jHHSD we conclude fp(t, X, S) = J l (t, X, S), i = 0, 1. 



(B.15) 
□ 



Appendix C. Numerical Scheme 

To solve for the indifference price p(t,S), we apply implicit finite-difference scheme to 
the system (|4.8p . Working with the log-price Z t = log St, we use an equally spaced mesh 
(l,Z), with grid points (ti,Zj), U = iS t , i = 0, . . . , N, S t = T/N and Zj = Z min + j5 z , 
j = 0, . . . , M, with 

= \l^t+(\° 2 ) 2 Sl (CI) 

Letting p l - = p(ti,Zj) denote the option price on the grid, and similarly for q % -, we use the 
standard second-order centered finite-difference approximation (see e.g. [25]) for all the 
first and second derivatives of p(t, S) and use one-sided approximations at the boundaries 
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Z = Z m i n and Z = Z max . To avoid the difficulty of the nonlinear term in the p-PDE, the 
following first order approximation is used 

e j P (U, Zj ) _ e ~/ P ) ^ e iP) +1 + ie iP) +1 _ p i+ly ( C-2 ) 

The resulting implicit tri-diagonal system for [p 1 -) is solved using the Thomas algorithm, 
see |24t Ch 3.5]. In the second half-step, we apply an explicit Euler scheme for the g-ODE 
(second line of (|4.8|) ) using the just-computed pj, 

which can be explicitly solved for q 1 -. 



Appendix D. General Illiquid Dynamics 

One could choose to work with ()2.5|) . the general case when \i\ ^ 0, o\ ^ 0, and L^O. 
In that case the redefined value functions U°(t, X, S) and U l (t, X, 5, 7r) = e~ lX R 1 ^, S, 7r) 
still allow separation of wealth, however since tt is non-constant during the liquidity shock, 
R 1 solves a two-dimensional PDE in (S, it) that in turn creates further 7r-dependence in 
the HJB equation satisfied by R°. With the extra state variable 7r, we no longer have a 
closed-form solution for the Merton problem (i.e. when h(S) = 0) either. In analogue to 
the indifference price definition in (14. 7p we can show that the indifference prices satisfy the 
HJB equations (cf. (Oj) ) 

Pt + ^ 2 S\ P ss - ipD + m>S PS + (D.l) 

2 7 7 i<o 

+ sup|/, 7r - ia V 7 - 7^PS7T - ^l^P e ^ e -^~L)S,(^L),)) = Q) 

% + ^OiS 2 {qss - iq 2 s) + ^cr 2 ir 2 (q n7T - jql) + a\-nS{q slT - jqsq^) - wfirSqs 

2 2 c i i 1 2 2 ^10 ^0 -7(0-0) , ^10 1 ^1 n 

2 7 fj 7 7^1 

with terminal condition: p(T,S) = q(T,S,n) = h(S). No analytic insights into (jD.ip seem 
possible. If we assume a static illiquid state = 0, 01 = 0), but L / 0, then the resulting 
q( L ) = q( L \t, S) is independent of tt (because ditt = 0), and the price equations (|D.1|) can 
be simplified to 

J L ) _1_ 1 ^ rf/Ji) ^J^N + „ ^ 4- m 1 

Pt + ~°o^ - 7(Ps J J + Mo-aPs + — - ~ r L) 

7 1 P 

+ supj/xovr - i0 o 2 vr 2 7 - ^SpfK - e 7P« e 7^ e -7<^(*,(l-^)} = ^ 

w _ ^^ e _ 7{P (,)_ g (,)) + ^10 _ 1 (^y = (D 2) 

7 ' 7 7 * 
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In <|H2]> the Merton terms F (L) = F^ L \t) and f[ L) = f[ L) (t) are free of vr as well, and 
are given implicitly via 



dF (L) 



dt 



v 01 F^ - bu P { W F W - ^ VifV - ^^(t)} = 0, 



+ ^(if - if ) = 0, F (L) (T) = if^(T) = 1. (D.3) 

We hope these derivations convince the reader that this more general route is fraught with 
tedious complexity and hence provide intuition for the use of our main model in Section [2j 



References 

1. D. Becherer, Rational hedging and valuation with utility-based preferences, Ph.D. thesis, Technical Uni- 
versity of Berlin, 2001. 

2. , Utility-indifference hedging and valuation via reaction-diffusion systems, Proceedings of the 

Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 460 (2004), 
no. 2041, 27-51. 

3. C. Blanchet-Scalliet, N. El Karoui, and L. Martellini, Dynamic asset pricing theory with uncertain 
time-horizon, Journal of Economic Dynamics and Control 29 (2005), no. 10, 1737-1764. 

4. F. Delbaen, P. Grandits, T. Rheinlander, D. Samperi, M. Schweizer, and C. Strieker, Exponential 
hedging and entropic penalties, Mathematical Finance 12 (2002), no. 2, 99-123. 

5. P. Diesinger, H. Kraft, and F. Seifried, Asset allocation and liquidity breakdowns: what if your broker 
does not answer the phone?, Finance and Stochastics 14 (2010), no. 3, 343-374. 

6. R.J. Elliott, L. Chan, and T.K. Siu, Option pricing and Esscher transform under regime switching, 
Annals of Finance 1 (2005), no. 4, 423-432. 

7. H. Follmer and M. Schweizer, Hedging of contingent claims under incomplete information, Applied 
Stochastic Analysis, Stochastics Monographs (M.H.A. Davis and R.J. Elliott, eds.), vol. 5, Gordon and 
Breach, London/New York, 1991, pp. 389-414. 

8. M. Frittelli, The minimal entropy martingale measure and the valuation problem in incomplete markets, 
Mathematical Finance 10 (2000), no. 1, 39-52. 

9. P. Gassiat, F. Gozzi, and H. Pham, Investment/ consumption problem in illiquid markets with regimes 
switching, 2011, Arxiv preprint arXiv:1107.4210. 

10. V. Henderson and D. Hobson, Utility indifference priemg-an overview, Indifference Pricing: Theory 
and Applications (R Carmona, ed.), Princeton University Press, 2009. 

11. R. Houssou and O. Besson, Indifference of defaultable bonds with stochastic intensity models, Arxiv 
preprint arXiv:1003.4118 (2010). 

12. A. Ilhan, M. Jonsson, and R. Sircar, Portfolio optimization with derivatives and indifference pricing, 
Indifference Pricing: Theory and Applications (R Carmona, ed.), Princeton University Press, 2009. 

13. I. Karatzas and S.E. Shreve, Brownian motion and stochastic calculus, vol. 113, Springer- Verlag, 1991. 

14. A. Kirilenko, A. P. Kyle, M. Samadi, and T. Tuzun, The flash crash: The impact 
of high frequency trading on an electronic market, Tech. report, Available on SSRN 
http: //papers . ssrn.com/sol3/papers . cfm?abstract-id=1686004 , 2011. 

15. T.S.T. Leung, A Markov-modulated stochastic control problem with optimal multiple stopping with ap- 
plication to finance, Decision and Control (CDC), 2010 49th IEEE Conference, IEEE, pp. 559-566. 

16. F.A. Longstaff, Financial claustrophobia: Asset pricing in illiquid markets, Tech. report, National 
Bureau of Economic Research, 2004. 

17. M. Ludkovski and H. Min, Illiquidity effects in optimal consumption-investment problems, Arxiv 
preprint arXiv:1004.1489 (2010). 

18. O. Menkens, Crash hedging strategies and worst-case scenario portfolio optimization, International 
Journal of Theoretical and Applied Finance 9 (2006), no. 4, 597-618. 

19. Y. Miyahara, Canonical martingale measures of incomplete assets markets, Probability theory and 
mathematical statistics (Tokyo, 1995) (1996), 343-352. 

20. R. Rouge and N. El Karoui, Pricing via utility maximization and entropy, Mathematical Finance 10 
(2000), no. 2, 259-276. 

21. E.S. Schwartz and C. Tebaldi, Illiquid assets and optimal portfolio choice, Tech. report, National Bureau 
of Economic Research, 2006. 



EUROPEAN OPTION PRICING WITH LIQUIDITY SHOCKS 



25 



22. M. Schweizer, Approximation pricing and the variance- optimal martingale measure, The Annals of 
Probability 24 (1996), no. 1, 206-236. 

23. T.K. Siu and H. Yang, Option pricing when the regime-switching risk is priced, Acta Mathematicae 
Applicatae Sinica (English Series) 25 (2009), no. 3, 369-388. 

24. J.C. Strikwerda, Finite difference schemes and partial differential equations, 2 ed., SIAM, Philadelphia, 
2004. 

25. P. Wilmott, S. Howison, and J. Dewynne, The mathematics of financial derivatives, Cambridge Uni- 
versity Press, Cambridge, 1995, A student introduction. 

26. T. Zhou, Indifference valuation of mortgage-hacked securities in the presence of prepayment risk, Math- 
ematical Finance 20 (2010), no. 3, 479-507. 

Department of Statistics and Applied Probability, University of California,, Santa Bar- 
bara, CA 93106-3110, USA. 

E-mail address: ludkovski@pstat.ucsb.edu, shen@pstat.ucsb.edu 



